Lattice calculation of low energy constants with Ginsparg- Wilson type fermions 
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We present a quenched lattice calculation of low energy constants using the chirally improved Dirac 
operator. Several lattice sizes at different lattice spacings are studied. We systematically compare 
various methods for computing these quantities, using pseudoscalar and axial vector correlators. 
We find consistent results for the different approaches, giving rise to f n = 96(2) (4) MeV, Jk = 



106(1)(8) MeV, f K /.U = l.H(l)(2), E 
m = 4.1(2.4) MeV and m s = 101(8) MeV. 

PACS numbers: 11.15.Ha, ll.10.Kk 



-(286(4)(31) MeV) J , the average light quark mass 



I. INTRODUCTION AND MOTIVATION 

Low energy theorems were derived quite early in the 
study of hadrons. When QCD was accepted as a prime 
candidate for the underlying quantum field theory it be- 
came clear that its (approximate) chiral flavor symmetry 
is spontaneously broken. One may use the correspond- 
ing ground state as a starting point for a systematic ex- 
pansion (Chiral Perturbation Theory, ChPT) taking into 
account the explicit symmetry violations due to small 
masses of the light quarks. In such an expansion phe- 
nomenological parameters, the so-called low energy con- 
stants, have to be determined independently. 

The leading order low energy parameters like the pion 
decay constant or the quark condensate are well-known 
from "classical" low energy theorems and experiments. 
It is a challenge, however, to find these parameters based 
exclusively on ab-initio calculations for QCD. Of course 
also QCD has its minimal set of input parameters (fixing 
the scale and the quark masses), but except for these, 
all other properties can, at least in principle, be de- 
rived. Due to the special nature as a strongly inter- 
acting quantum field theory only non-perturbative cal- 
culations lead to that numbers. The lattice formulation 
appears to be the prime approach for that aim, using a 
non-perturbative computer evaluation of the QCD path 
integral. 

For a long time the lattice approach has been plagued 
by the problem of incorporating chiral symmetry in the 
fcrmionic action. It took almost two decades from the 
first numerical results for lattice gauge theory before it 
was understood how chiral symmetry may be formulated 
on the lattice. The Ginsparg- Wilson condition (GWC) 
PJ , characterizes the class of Dirac operators allowing for 
the lattice analog of chiral symmetry H. The Dirac op- 
erators obeying the GWC exactly jj, US @ or approx- 
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llll | allow one to approach smaller 
quark masses than with the more traditional Wilson 
fermions. 

However, the GW-fermions need additional effort 
for their numerical implementation. Overlap fermions, 
which obey the GWC exactly, are almost two orders of 
magnitude more expensive in terms of computer power 
than Wilson fermions. The approximate GW-fermions 
are cheaper, but still one order of magnitude more ex- 
pensive. For that reason there are only few first attempts 
to use these Dirac operators for a simulation of full QCD 
with dynamical quarks. There are, however, several re- 
sults for the quenched case. 

In ^3 hadron masses in the quenched case have been 
studied for two types of Dirac operators that obey the 
Ginsparg- Wilson relation to a good approximation: the 
fixed point operator and the chirally improved operator. 
Due to their good chiral behavior, these fermions provide 
a suitable framework for the calculation of low energy 
constants. Since such a calculation involves quantities 
that are renormalized, one also has to determine renor- 
malization constants for the connection with a contin- 
uum scheme like MS. Whereas in the case of the overlap 
action there are exact symmetries relating the renormal- 
ization constants of the scalar with the pseudoscalar and 
the vector with the axial vector sectors, for fermions with 
only approximate GW-symmetry such relations have to 
be checked, too. 

For the overlap action there have been several deter- 
minations of low energy parameters 0, Q, 0, 0, 0, 
13 IH HI EJ IH EI |1 |2|. For the chirally im- 
proved Dirac operator Dqi only preliminary results have 
been published 0,E(j. The reason was, that the renor- 
malization constants were not available. Meanwhile the 
necessary constants for quark bilinears have been deter- 
mined for Dq\ in |27j. This now allows us to compute 
some of the basic low energy parameters in the quenched 
case. For this purpose we study quenched QCD at vari- 
ous values of the quark masses and determine results for 
m u = nrid ^ m s in the ud and us meson sector for several 
lattice sizes and lattice spacings, down to a pion mass of 
330 MeV. Our gauge action is the Liischer-Weisz action 
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In Sect. ^ we fix our notation and recapitulate the ba- 
sic relations like PCAC, the axial Ward identity and the 
Geil'Mann-Oakes-Renner (GMOR) relation. We then 
discuss the observables used in our analysis. In Sect. IIIII 
we discuss the details of the lattice simulation and present 
the results for the masses and low energy parameters (like 
the meson and quark masses, meson decay constants, chi- 
ral condensate) in Sect. lIVl We summarize and conclude 
in Sect. E| 



II. LOW ENERGY CONSTANTS AND THEIR 
RELATIONS 

All of the relations are given for Euclidean space-time. 
We first introduce our notation and then briefly summa- 
rize the three key relations which we need in our extrac- 
tion of low energy constants from lattice simulations. A 
collection of the relevant vacuum expectation values and 
ratios thereof ends this section. 



A. Basic definitions and relations 

For flavor symmetry group SU(2) we define the un- 
renormalized (lattice) isovector operators for the pseu- 
doscalar, vector and axial vector sector through (all fields 
are taken at the same space-time point) 

P a = ^75T a V, (1) 
V = ~?7„t<^, (2) 
A; = \ jn^r a ^. (3) 

The fermion fields are flavor doublets ip = (u, d) and 
r a ,a = 1,2,3 are the Pauli-matrices. Charged quark 
bilinear operators are defined through, e.g., 

P+ = (r 1 -it 2 ) 4; = d 75 u . (4) 

Quantities renormalized according to a continuum 
renormalization scheme are denoted with a superscript 
(r). As our reference scheme we use the MS-scheme at 
a scale of /i = 2 GeV. The renormalized quantities are 
related to their lattice counterparts via renormalization 
constants, 



p( r ) a 


= Z P P a , 


(5) 


Tr(r)a 


= z v v; , 


(6) 


A(r)a 
A< 


= Z A A«, 


(7) 




= Z m m , 


(8) 




= Z S Y> with Z s = 1/Z m . 


(9) 



Here denotes the renormalized quark mass and EM 
the renormalized condensate, where the bare condensate 
reads 

£ = -(uu + dd) . (10) 

The renormalization constants relate different schemes. 
For exactly chirally symmetric actions we would have 
Zg — Zp and Zy — Z A . The chirally improved Dirac 
action used here is only approximately chirally symmet- 
ric and we have com put ed the renormalization factors 
(in the chiral limit) in 27], utilizing the non-perturbative 
methods suggested in [23, H3| • 

We now discuss the basic relations in terms of renor- 
malized quantities (measurable experimentally and de- 
fined in the MS-scheme). The axial vector current is re- 
lated to the vector current by commutation relations in 
current algebra. For conserved vector currents the nor- 
malization is therefore fixed for both. The axial vector 
operator couples to the weak interaction currents. Its re- 
lation to the physical (renormalized) isovector pion field 
defines the pion decay constant, i.e., 

d^ a = Ml U4> [r)a ■ (ii) 

(There are also other conventions differing by, e.g., a fac- 
tor of Our definition corresponds to an experimental 
value of 92.4(3) MeV The (renormalized) pion field 

obeys (x = (a?, t)) 

(O|0(x)« a \tt\p = 0)) = ^e"^' . (12) 

Various equivalent expressions for expectation values 
may be derived, in particular the so-called partially con- 
served axial vector current (PCAC) relation: 

a M <0| 4'>(x) \7r b (p - 0)) = S ab Ml U e~ M ~ 1 . (13) 

Global symmetries of an action at the classical level 
lead to conserved Noether currents. In a quantum field 
theory global symmetries of the action are expected to 
manifest themselves by an analog to the conservation of 
Noether currents. To derive these relations one considers 
infinitesimal symmetry transformations of the fermion 
fields in the path integral 

(0|C|0) = i J dU dip dip 0[U,4>,ip] e - s[u ^^' ] , (14) 

where O denotes an arbitrary operator. Some classical 
symmetries may be spontaneously broken or broken due 
to anomalies resulting from the functional integration. If 
the integration measure is invariant under the transfor- 
mation this lead to (Ward-) identities of the form 

(0\8O\0) - (0\OSS\0) = . (15) 

Thus, the quantum analog to the classical conservation 
law has similar (or equal, O = 1) form, but is an operator 
identity. 
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For the QCD Lagrangian (renormalized, e.g., in the 
MS-scheme) one exploits the invariance properties under 
symmetry transformations to derive Ward identities. Al- 
though obtained first on the classical level and on-shell, 
the relations hold under quantization and one ends up 
with local operator identities for the full quantized the- 
ory. (For the singlet axial vector there is an additional 
anomaly contribution from the integration measure.) 

In the quark sector chiral symmetry is broken explicitly 
by the quark mass matrix M. and one finds 



\ (V4r a ,.M]V>) (r) . (16) 
~ (MT a ,M}^f ] . (17) 



If the quark masses are degenerate, the vector current is 
conserved and the axial vector current obeys the axial 
Ward identity (AWI), 

9 M 4 r)a = 2m(r)p(r)a - (18) 
Combining this with relation gives 

2m' r) F (r »" = / I Jl4 2 f • (19) 

The correlation function of the normalized pion field 
(cf. (T2J) reads 

(0\4><r> a (j?= O,i)0 (r)b (O)|O) = -L-8a b e- M - t . (20) 

The asymptotic behavior of the pseudoscalar field corre- 
lator 

(0|P«>- = 0,i)P« b (0)|0) = (G«) 2 -L- 



Sabe 



-M„ t 



. ( 21 ) 

is dominated by the pion state as well. We have intro- 
duced the relative factor , which, from (|19|l , is 



\(0\P^\n(0))\=G^ = 



(22) 



2toM 

Another Ward identity may be derived taking O = P, 
leading to 



E (r) 



Ml 



|(o|p«K(o))| 2 , 



(23) 



where we assume degenerate light quark masses for sim- 
plicity. This, combined with the definition of the pion 
decay constant and 119(1 , leads to the usual form of the 
GMOR relation: 



(24) 



The exploitation of the underlying principles has then led 
to the development of chiral perturbation theory |32l l33| . 
In that context a systematic expansion of many observ- 
ables in terms of low energy constants has been derived. 
These constants, however, have to be determined ci- 
ther from experiment or from basic principles, i.e., non- 
perturbative solution of the underlying field theory QCD. 
The lattice formulation of QCD allows such a determi- 
nation. 



B. Lattice relations 

Let us now express the renormalized quantities 
through their lattice counterparts, using the relations 
|[SJ - 10. Since we need to avoid the evaluation of dis- 
connected diagrams we restrict ourselves to the isovector 
charged bilinears 



-^(t 1 -ir 2 )V = dTu , 



(25) 



where Y is 75, 7^ or 7^75. Consequently from now on we 
drop the flavor superscript a. 

In the simulations we always keep the quark sources 
(and thus the operator source) at a fixed time t = and 
compute propagators to all other lattice points. We sum 
over the spatial volume of the sink time slice in order to 
project to zero spatial momentum p = 0. 

Taking into account the renormalization factors, the 
following correlators of lattice operators and ratios have 
been studied (ji = 4 refers to the Euclidean time direc- 
tion, the symbol ~ denotes the asymptotic behavior for 
large t): 



Z 2 P (P(p = 0,t)P(0)) 



-M„ t 



Z A Z P (A 4 (p = 0,t)P(0)} 



f 2 M 3 

4( TO M)2 

M„|S( r >| 
f 2 M 2 

| S M| e -A/^ 



e 

-Af* t 

■M, t 



,(26) 



(27) 



Z A Z P (d t A 4 (p = 6, t)P(O)) 



f 2 M 3 



-M„ t 



7 2 



z 2 A {d t Mp 



A 4 (p = 6,t) A 4 (0)) 
6,t) At(0)) 



2 m( r ) 

|E< r >| Af^e - -"-* (28) 
M K fie- M **, (29) 
M 2 f 2 e~ M ^, (30) 



Z A Z P ^J(A 4 (p = 6,t)A 4 (0)) (P(p = 5,t)P(0)> 



-M„ t 



(31) 



The pseudoscalar masses thus may be derived from the 
exponential decay and the other low energy parame- 
ters from its coefficient. We cannot use correlators like 
((<9 t A4(i))(<9t Ai(0))) since the source is fixed to the time 
slice t = and thus we cannot construct the lattice 
derivative there. 
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The asymptotic behavior cancels in the following ra- 
tios: 



Z A (d t A 4 (p = d,t)X(0)) 

Zp (p(f=6,t)x(o)) 



r (r) 



= Z m 2m = 2m (r) , (32) 
where X may be P or A 4 . Further useful ratios are 



Z A (A 4 (p = 0,t)P(0)) 
Z P (P(p=0,t)P(0)) 



Z\ (A 4 (p = 0,t)A 4 (0)) 

z% (p(f=6,t)P(o)) 



Z P (d t A 4 {p = 6,t)P{0)) 
Za (A 4 (p = 0,t)A 4 (0)) 



2m 



2to« M^fl 



r (r) 

2m' r) 

r (r) 



(33) 



Ml 



(34) 
(35) 



TABLE I: Parameters of the simulation. Where the strange 
quark mass is given (in parenthesis) , we also determined prop- 
agators for strange hadrons. Type denotes the type of the 
quark source/sink and #cf. the number of configurations en- 
tering the analysis. 



L 3 


x T 


P 


a [fin] 


alGeV" 1 ] 


#cf. 


Type 


a m 


(am s ) 


8 3 


x 24 


7.90 


0.148 


0.750 


200 


p,n 


0.02 


- 0.20 


12 3 


x 24 


7.90 


0.148 


0.750 


100 


p,n 


0.02 


-0.20 


12 3 


x 24 


8.35 


0.102 


0.517 


100 


p,n 


0.02 


-0.20 


16 3 


x 32 


7.90 


0.148 


0.750 


99 


p,n 


0.02 


-0.20 


16 3 


x 32 


8.35 


0.102 


0.517 


100 


p,n 


0.02 


-0.20 


16 3 


x 32 


8.70 


0.078 


0.395 


100 


p,n 


0.02 


-0.20 


16 3 


x 32 


7.90 


0.148 


0.750 


100 


p, n, w 


0.02 


-0.20 
















(0.08,0.10) 


20 3 


x 32 


8.15 


0.119 


0.605 


100 


p, n, w 


0.017 


-0.16 
(0.06) 



TABLE II: Renormalization constants taken from |2jj. The 
values for /3 = 8.15 have been obtained by interpolation. 



Z s 



Zt 



Za 



zSub 



7.90 1.1309(9) 0.9586(2) 0.9944(3) 1.0087(4) 1.0281(5) 

8.15* 1.081(1) 0.967(1) 1.014(1) 1.011(1) 1.012(3) 

8.35 1.039(1) 0.973(1) 1.028(2) 1.012(1) 0.987(4) 

8.70 0.959(2) 0.979(1) 1.049(1) 1.0095(7) 0.915(1) 



(d t A 4 (p = 6,t)p(o)) 2 M2 

(A 4 (p=0,t)A i (Q))(P(p = Q,t)P(0)) 

Ratios involving the lattice derivative dtA 4 depend on 
the way the derivative is taken. Details will be discussed 
m Subsection HITT1 

Due to lattice periodicity and the parity properties of 
the meson propagators, the exponential term exp (— Mt) 
is accompanied by another term from the propagation 
backwards in time, exp (— M (T — i)). Depending on the 
correlation function we therefore observe cosh- (for (PP) 
and (A 4 A 4 ) correlators) or sinh-behavior (for (A 4 P) and 
((d t A 4 )A 4 ) correlators). 

III. TECHNICALITIES 
A. Setup 

In the fermion action we use the chirally improved 
Dirac operator Dqi 0, [nj. It is based on a system- 
atic expansion of the lattice Dirac operator taking into 
account the whole Clifford algebra and terms coupling 
fermions within a certain range of neighbors on the lat- 
tice. The expanded Dirac operator is inserted in the 
GWC which then leads to a set of algebraic equations 
for the expansion coefficients. The equations include a 
normalization condition, which depends on the lattice 
spacing and thus on the gauge action. The solution of 
the system of algebraic equations gives the coefficients 
defining Dei- In recent applications as well as here, 
we use a set of 19 independent terms in the expansion in 



the action. Taking into account the lattice symmetries 
this then corresponds to several hundred coupling terms 
at each site of the lattice. Finally, in the definition of 
Dei also one step of HYP-smearing of the gauge con- 
figuration is included (34|. For the gauge fields we use 
the Luscher-Weisz action [2^. The lattice spacings have 
been determined using the Sommer parameter [3Jj|. We 
summarize the simulation parameters in Table [I] 

The strange quark mass parameter was set using the 
kaon mass (chiral extrapolation with non degenerate 
quark masses). This gives a value |36| of am s — 0.089(2) 
for /3 = 7.90 and am s = 0.061(1) for (3 = 8.15. For 
P = 8.15 we have quark propagators at am = 0.06 and 
we used those for the strange quark. For (3 = 7.90 we 
have quark propagators for am = 0.08 and am = 0.10. 
We interpolate the hadron results from those two values 
to the strange quark mass at am = 0.089. All results 
for the if-meson are based on quark propagators with 
non-degenerate quark masses. In the discussion the cor- 
responding meson fields, however, will still be called P 
and A in order to simplify our notation. 

The necessary renormalization constants for fcrmionic 
bilinears were computed in |27| and we list them in Ta- 
ble [H] for completeness. All Z are in the chiral limit and 
for the MS-scheme at a scale of 2 GeV. The values at 
[3 = 8.15 were obtained by interpolation. 
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B. Sources and normalization 

In order to improve signals we use Jacobi-smeared 
[3?],|38| sources at t = for the quarks. For most runs we 
have one smearing width (denoted as n(arrow), cf. Ta- 
ble!]}. For some data sets, that have also been used for an 
analysis of the excited hadron masses [39l | , we have two 
different widths (n and w(ide)). In addition we also use 
point-like sinks (p(oint)) for sake of normalizing our re- 
sults as discussed below. The smearing widths are chosen 
such that the effective size of a given source is approxi- 
mately the same for all lattice spacings. The parameters 
(smearing steps and coupling) thus depend on the gauge 
coupling [39|. 

For computing the coefficients of the exponential decay 
necessary for determining the pion decay constant and 
the condensate, we have to normalize the hadron sources 
to point-like bilincars. Our normalization procedure is 
set up as follows: Let us denote by X SlS2 a mesonic op- 
erator built out of an anti-quark of smearing type si and 
a quark of smearing type S2 with Sj = n, w or p. Guided 
by the arguments discussed below we find that we may 
safely assume factorization of the normalization factors, 



X Sl s 2 — cf 2 x pp 



For ratios like 



(X SlS2 X s 
(X SlS3 X. 



(37) 



(38) 



wi 



we find excellent plateaus for t, t' G [t a , h] = [4, T — 4] 
and violations of the factorization hypothesis (|37l) smaller 
than 1.2 % for the narrow and wide results at j3 = 7.90, 
and less than 2.2 % for the results at j3 = 8.15. In this 
way we obtain numbers for (Cjf ) 2 , (C^f) 2 and C^C* 
and from them find the values for , . The coeffi- 
cients are calculated for all operators X, and separately 
for each gauge coupling, volume and quark mass. For the 
large physical volumes we find very little 0(1%) depen- 
dence of the on the quark mass which increases for 
the smallest volumes to a variation of 0(5 — 10%) over 
the range of quark masses we consider. The relative sta- 
tistical errors from the plateau fits are less than 0.0005 
and were neglected in the further analysis. 

As a consistency check we compared the final re- 
sults for all possible smearing combinations and checked 
whether the resulting masses and in particular the pref- 
actors of the exponential decay are in agreement when 
normalized appropriately. This was the case and in the 
subsequent presentation we therefore do not discuss this 
issue any more. All results presented here are based on 
the results with narrow quark sources. 



C. Fits, error estimates and numerical derivatives 

Fitting procedure: At larger t, where excited state 
contributions become negligible, the propagators have a 



cosh or sinh functional behavior. Therefore we may ex- 
tract the prefactor D and the meson mass M performing 
a correlated least-squares fit of the correlation function 
C{t) to 



- f e -Mt ± e -M (T-t) 



D(M)f(M,t) with /(M,t)=(. 
by minimizing 

X 2 = ^2[c(t)-D(M)f(M,t)]Cov-\t,t') 



(39) 



[C(f ) - D(M) f(M, 



(40) 



Here Cov denotes the covariance matrix for the correla- 
tion function entries and t, t' G [t a ,£f,]. The minimiza- 
tion may be simplified by observing, that for given M 
the minimum of \ 2 is obtained for 



D{M) 



E M , C(t) Cov' 1 (t, t<) /(M, f) 
Et,t' /( M > t) Cov' 1 (t,t')f(M,t>) 



(41) 



Using this relation one performs the one-dimensional 
minimization of l|40|) with regard to M. 

Ratios of correlation functions are fitted to constants 
in t e [t a ,h]. 

Error estimation: The correlator values are arranged 
in overlapping blocks according to the standard jackknife 
algorithm. Each block consists of typically 95% of all 
hadron propagators for a given set of lattice parameters. 

For each such jackknife block we then determine the 
values of the propagator and the requested ratios at var- 
ious i, as well as the covariance matrix for the propaga- 
tor and the variance for the ratios. These are then fit- 
ted as discussed, i.e., either to asymptotic cosh- or sinh- 
behavior or to a constant. For this fit the relative weights 
(as defined from the covariance and variance) are impor- 
tant. 

We also need variances for the correlators involving 
derivatives with regard to t (discussed below), as they 
also enter some of the ratios. Since we need these for each 
jackknife set we estimate the variance of the derivatives 
by performing another jackknife analysis within the given 
set. 

The fits (and the extrapolations to the chiral limit) are 
then repeated for all jackknife blocks and the variation of 
the results for coefficients, mass values and ratios leads 
to the estimate for the corresponding errors. 
Numerical derivatives: For some of the ratios of cor- 
relators we need derivatives of the correlator with regard 
to time t. Numerical derivatives are always based on as- 
sumptions on the interpolating function. Usual simple 2- 
or 3-point formulas assume polynomials as interpolating 
functions. We can do better by utilizing the information 
on the expected cosh- and sinh-dependence. In fact, we 
may use these functions for local 3-point interpolation 
and get the derivative therefrom. 
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A 2-point derivative based only on function values yt 
and yt-i is not suitable since it provides values at half- 
integer t. We therefore use a local 3-point fit to yt-i, 
yt and yt+i to the functional form l)39|) (depending on 
whether the correlator is symmetric or anti-symmetric in 
t), where the parameters D and M now depend on the 
actual value of t. We then reconstruct the derivative as 



d t f(M, t) = M 



-Mt 



±e 



-M(T-i) 



(42) 



Then the desired ratios can be computed in a straight- 
forward manner and the plateau values obtained by fits. 

When analyzing not ratios but the correlators of type 
((dtAi)X) directly, we perform a global fit to (A4X) ac- 
cording to H39|) and take the analytic derivative. 



IV. RESULTS 

Throughout the analysis we refer to the quantities as 
defined in Sect. Ill 61 Where we give physical, renormal- 
ized (MS-scheme) values we always use the renormaliza- 
tion factors as obtained in the chiral limit, even in plots 
for non- vanishing bare quark masses. The rcnormaliza- 
tion factors Z relate to the continuum MS-scheme at a 
scale of fj, = 2 GeV and are a-dependent as given in Ta- 
ble nn 



A. Meson Masses 

In Fig.m we present results for the pseudoscalar 
masses as determined from the exponential fit of the A4- 
correlator (|29|l for the time components of the axial vec- 
tor. The kaon propagators are determined from the same 
interpolators but for non-degenerate quark masses. The 
translation of results in lattice units to physical mass val- 
ues was done using the values for the lattice constant a 
given in Table HJ It is interesting to note that the linear 
extrapolation of our pion mass results implies a small 
residual quark mass. We remark that when using the 
pseudoscalar correlator at f3 = 7.90 and an extrapolation 
formula inspired by q uenched chiral perturbation theory 
(QChPT) jHElEJ no such effect was observed H^. 

The value of the strange quark mass has been cho- 
sen such that the kaon mass agrees with experiment at 
the linear extrapolation in the light quark mass to the 
physical point. For the 16 3 x 32 lattice at (3 = 7.90 
we interpolate M\ linearly between the two values of 
the strange quark mass neighboring the physical point at 
am s = 0.089, cf. Fig. El 

Let us now discuss possible sources of systematic er- 
rors. 

Finite size dependence: In Fig.|3 we compare 
a 2 G^ ] /Z P determined from ® for am > 0.01 for 
three volumes but at the same lattice scale. We find 
a strong volume dependence for the smallest lattice at 
bare quark masses below 0.06/a. This is to be expected 
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FIG. 1: M% and M K determined from the asymptotic behav- 
ior of the ^-correlator 129H . 
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FIG. 2: M K vs. M n . For the 16 3 x 32 lattice we have in- 
terpolated between two values of the strange quark mass; for 
the 20 3 x 32 lattice am s — 0.06 corresponds already to the 
physical point. The curves represent the linear fit in Fig.0 



since for that quark masses the inverse pion mass be- 
comes larger than half the spatial lattice size. For the 
larger lattices finite volume effects are significant only 
for am < 0.02. We will therefore discuss only results 
for the large lattices and am > 0.02 (0.017 for the 
20 3 x 32 lattice). The lattices largest in physical units 
and of similar spatial extent 0(2.4 fm) are those of size 
16 3 x 32 at (3 = 7.90 (a = 0.148 fm) and 20 3 x 32 at 
(3 = 8.15 (a = 0.119 fm). 

Topological finite size effects: Exact Ginsparg- 
Wilson Dirac operators have exact zero modes. In 
quenched simulations these are not suppressed by the 
fermionic determinant. However, their effects for, e.g., 
the pion correlators |43|. are hard to detect unless one 
approaches very small pion masses 0(250 MeV) (SJ- For 
approximate GW-operators like the one studied here the 
main problem is the (scarce) occurrence of slightly mis- 
placed zero modes, real eigenvalues of Dq\ below zero. 
These lead to contributions in the quark propagator that 
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FIG. 3: a 2 Gi r) /Z P from ^ for = 7.90 and lattice sizes 
8 3 x 24, 12 3 x 24, 16 3 x 32. The finite size significantly affects 
the small volume data (8 3 x 24) for am < 0.06. For the other 
volumes finite size effects become important for am < 0.02. 
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FIG. 4: Ratio of the pion masses obtained from the pseu- 
doscalar propagator (PP) and the axial propagator {A4A4) 
for the same physical volumes but two different lattice spac- 
ings. 



diverge already at some positive, albeit small, mass. The 
smaller the quark mass parameter is, the stronger one 
might notice such distortions often referred to as topo- 
logical finite size effects. 

There have been various suggestions to deal with the 
problem of zero modes in the pseudoscalar propagator 
ElllllliEli among them the proposal to study a com- 
bination of the iso-non-singlet pseudoscalar and scalar 
correlators such that the zero modes cancel. Indeed, zero 
modes contribute differently to different propagators. For 
exactly chiral Dirac operators for the ^-correlator we 
expect a contribution of 0(1/ m), whereas for the P- 
correlator we expect C(l/m 2 ). 

In Fig. 21 we show the ratio of pseudoscalar masses 
obtained from the pseudoscalar and axial propagator 
M„ ( ( P P ) ) / M„ ( ( A A A A ) ) , comparing results at two dif- 
ferent lattice spacings for the same physical volumes. 
There are indications for zero mode contributions being 
stronger for the P-correlator on the coarser lattice. For 
the final pseudoscalar mass values we use the less affected 
results from the ^-correlators. 

Chiral logs: Since quenched QCD does not include 
the full fermion dynamics, corrections to the chiral ex- 
pansion introduce new terms, including singular ones 
(quenched chiral logarithms). On the other hand, it has 
been hard to clearly identify these in numerical calcula- 
tions 

El 13 13 In particular, one needs to approach 
parameter regions, where the pion mass is well below 300 
MeV. The situation then is furthermore obscured by the 
role of zero modes (depending on the fermion action), 
poor statistics and variation in the fit range. 

In quenched chiral perturbation theory (QChPT) [U 
I4TL |43| one expects for the pseudoscalar mass non- 
analytic behavior in the quark mass parameter 0, OH 1 
which arises from hairpin diagrams and rj would-be- 
Goldstone bosons 

(dMpfoctam) 1 /^), (43) 



where values for 5 ranging from 0.19 to 0.23 have been 
quoted in recent literature 0, 0, 0] . In the range of 
values studied here we cannot disentangle reliably such 
effects from the possible spurious zero mode contribu- 
tions and consequently do not attempt to determine 8. 



B. Quark masses 

In Eq. (|32|l one utilizes the axial Ward identity to ob- 
tain the so-called AWI-mass toawi! the asymptotic expo- 
nential behavior of numerator and denominator cancels 
and thus, including the renormalization factors, this al- 
lows one to obtain the renormalized quark mass. In Fig.0 
we plot the results for two choices of the pseudoscalar in- 
terpolator P, showing the results both in lattice and in 
physical units. 

For the second choice, X = A4, the correlators are 
of sinh-type and therefore the ratio becomes numerically 
unstable near the symmetry point in t. We thus cannot 
reliably (with small errors) determine the plateau val- 
ues for this choice if the quark mass parameter is small. 
However, at all masses where we can compare the two 
results they are in excellent agreement as may be seen in 
the figure. The subsequent analysis is based on the more 
stable choice X = P. 

In Fig. we compare the values, again obtained utiliz- 
ing ll-'il't with X — P. The are derived for two different 
lattice spacings and given in physical units. The lines 
correspond to linear fits m^ r ' oc M% , enforcing the simul- 
taneous chiral limit of both observables. The data do 
not show deviation from that behavior, although loga- 
rithmic corrections are expected due to quenching. We 
observe that the linear extrapolations are in good agree- 
ment with the Particle Data Group [3lJ average for the 
light quark masses at the physical pion mass. Averag- 
ing the extrapolated values obtained at the physical pion 
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FIG. 5: The renormalized quark mass arrv r ' (in the MS- 
scheme) vs. the bare mass parameter am as determined from 
l|32l (16 3 x 32, p = 7.90) using X = P (squares) and X = A 4 
(circles). For X = A4 values at quark masses below am — 
0.06 have been omitted due to unstable plateau fits. 
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FIG. 6: Plot of m (r) (obtained from vs. M% for two 

different lattice spacings but similar physical spatial lattice 
volume. The linear extrapolations are in good agreement with 
the expected light current quark masses. 



mass we obtain an average light quark mass of 



5 K 



4.1(2.4) MeV 



(44) 



in the MS-scheme. The error is essentially due to the 
residual quark mass (see Fig. s El and [I}. The precise de- 
termination of this residual quark mass is obscured by 
the possible contribution of a quenched chiral log (cf. the 
discussion in 01 )• This effect is stronger for the smaller 
physical volumes and thus we refrain from a determina- 
tion of m( r ) for the lattices at /3 = 8.35 and 8.7. This 
prohibits a continuum extrapolation linear in a and we 
therefore quote the average of the two values for the two 
lattices with spatial extension 2.4 fm. 

As addressed in Sect. IIIII the bare mass parameter for 
the strange quark has been fixed such that the kaon mass 
assumes its physical value when the data are extrapo- 
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using X — P. The slope of 
m^ for the pion data provides Z m . 



TABLE III: Renormalization constant Zs from |27l| (see Table 
ITTt compared to the values of Z m as derived from the slope of 
the renormalized quark mass. 



P a [fin] o[GeV _1 ] Z s 



1/Zs 



ZmZs 



7.90 0.148 0.750 
8.15 0.119 0.605 



1.1309(9) 0.8842(7) 0.891(4) 1.007(5) 
1.081(1) 0.9250(9) 0.916(5) 0.991(6) 



lated to the physical pion mass |3j| . Fig. compares 
the results for the corresponding renormalized masses 
(in lattice units) derived from (|32[1 for kaon and pion 
correlation functions. For the light quark masses (the 
pion propagator) the results confirm those of Fig. El now 
for two different lattice spacings superimposed. For the 
strange quark the fitted lines are close to being paral- 
lel with negative intercepts, in good agreement with the 
negative bare strange quark mass parameters. This again 
demonstrates that our Dirac operator has very small ad- 
ditive mass renormalization. 

For the pion data, the slope of provides a value 
for the light quark mass renormalization factor Z m . In 
Table IIIII we compare these numbers with the values of 
XjZg from |27| and find very good agreement. From 
the PCAC-relation and l)32|l we expect for the kaon data 
am ( r ) ~ z m a (m s + m)/2 which is indeed observed. 

The quark mass data shown in Fig. [SI now in physical 
units in the MS-scheme, again give a very consistent pic- 
ture. Note that now the abscissa gives the corresponding 
pseudoscalar mass, i.e., that of the pion or the kaon, for 
the corresponding values of m^ r \ The error bands show 
that for given lattice spacing the numbers for kaon and 
pion are on top of each other, although these states have 
different quark content. (This justifies the often used 
strategy to obtain kaon results from pion propagators 
using degenerate mass quarks.) 

Like for the light quark masses in Fig. [HJ also the result 
for the renormalized strange quark mass is in excellent 
agreement with the world average [3l| . Averaging the 
two values computed at the physical kaon mass for the 
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FIG. 8: m (r) vs. Mp S from ^ using X = P. Here M PS FIG. 9: 



denotes M n or Mj-, respectively. All masses are given in 
physical units in the MS-scheme. 



two lattice spacings we obtain for the strange quark mass 



|E (r) | vs. am for 16 3 x 32, /3 



7.90 and 

20^ x 32, (3 = 8.15 comparing different way of extraction; 
circles: GMOR-slope Eq. g]J, boxes: {A 4 P} Eq. dia- 
monds: (A^A^^PP) Eq. (13 IB . We also show linear extrapo- 
lations to the chiral limit. 



- (m< r) + ™ (r) ) ~ 52(3) MeV (45) 

in the MS-scheme. With (|44|1 for the light quarks this 
then gives mi^ = 101(8) MeV. Possible finite size effects 
and other systematic effects like chiral extrapolation and 
quenching have not been accounted for. The given error 
takes into account the standard error and includes the 
derivations due to the dependence on the lattice spac- 
ing. As discussed above for the light quarks, we have 
only values at two lattice spacings available and there- 
fore cannot perform a sensible continuum extrapolation. 
Our numbers for the average are in good agreement with 
determinations from the overlap action in |l4l Il7| , and 
to a less extent also with I25L 



C. Condensate 

For exact GW-fermions the chiral condensate may be 
obtained from the trace of the inverse Dirac operator. In 
our case the subtraction constant is not known to suffi- 
cient precision and therefore this approa ch, successfully 
applied for the overlap operator |2lL |22| , does not work 
|47|. Another method uses the density distribution of 
the small complex eigenvalues; this would be a promising 
method in our case but it requires the costly determina- 
tion of the low lying eigenvalue spectrum for all configu- 
rations which we did not calculate. 

Instead we have computed the renormalized conden- 
sate from the relations H24fl . (|27|l and (|31H which all con- 
tain £M. The first of these is the GMOR relation, the 
other two are determinations directly from the coeffi- 
cients of propagators and implicitly related to GMOR 
as well. In Figs. ED and ^| we show our results for E^ r ' 
in lattice units as a function of the bare quark mass. 



I 

i 




FIG. 10: |E (r) | vs. Ml as in Fig.^l but in physical units. 

We find excellent agreement for all three determinations. 
The values are consistent within the error bars. 

The dependence on the bare quark mass is compatible 
with the leading linear chiral behavior. Note, that we are 
not at small enough quark masses to be in the so-called 
e-regime but are in the p-regime |48| (for recent studies 
in the e-regime cf. [H |Hl Eg). We show the linear fit 
(with 1 s.d. error band) to the data from all three types 
of determination but omit the points with smallest mass 
value in the fit (cf. our discussion of possible finite size 
effects). 

So far we have restricted ourselves to the discussion of 
the largest physical volume results for only two lattice 
spacings, corresponding to — 7.90 and 8.15. In order 
to analyze the scaling behavior in the continuum limit we 
show in Fig.^Jthe values for all lattice spacings (always 
for the largest volume) studied. Different physical lattice 
volumes are combined in this figure. As may be seen from 
the upper part of the plot, the leading a 3 -dependence is 
suggestive. 

In the lower part of the figure we exhibit the renor- 
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FIG. 11: a 3 |E (r) | and |S (r) | vs. a. In the legends we indi- 
cate the quantities used in the derivation of the condensate 
value. In the bottom plot we show a constant extrapolation 
(horizontal error bands) and compare it to a linear one. 



malized condensate by dividing out the expected lead- 
ing 0(a 3 ) scale dependence. Since the Dirac action used 
obeys the GWC not exactly we cannot exclude linear cor- 
rections to perfect scaling. In these were found to be 
small for the hadron masses. For the condensate some of 
these corrections are already taken into account by the 
(independently determined) Z-factors in the determining 
equations. 

For estimating the systematic error of the contin- 
uum extrapolation we therefore perform both, a con- 
stant fit (no scaling corrections) and a linear fit to the 
a-dependence. We do this for each type of derivation. 
The average of the resulting values for the constant fit 
is = (286(4) MeV) 3 ; the linear fit leads to a larger 

value |E( r )| w (318(25) MeV) 3 . For the final value we 
give the result of the constant fit but quote the differ- 
ence to the linear fit result as systematical error. Thus 
we find 



= (286(4) (32) MeV) E 



(46) 



This value is slightly larger, although still in the error 
limits of a determination from the overlap action llal 
and larger than the corresponding results in Ref. P"t1.|24|. 
Our numbers are in good agreement with calculations 
in the e-regime 0, lljj . They also agree with a recent 
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FIG. 12: The dimensionless decay constants a f n and a fx, 
determined with Eq. 12911 . vs. am. 
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continuum extrapolation of the condensate based on the 
spectral distribution of the overlap operator j2^| • 



D. Pion decay constant 

The pseudoscalar decay constants have been extracted 
from the asymptotic behavior of the pseudoscalar cor- 
relation functions according to (|29|l for pion and kaon, 
respectively. (See also (Boj for another method to deter- 
mine /„.) 

In Fig. E| we compare these results. When plotting 
them as functions of the respective pseudoscalar masses 
in Fig. 1131 the data for pion and kaon essentially overlap 
each other and exhibit a universal functional behavior. 
We also show the error band of a quadratic extrapolation 
to the chiral limit. 

For full QCD the chiral expansion of the pion decay 
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constant should behave like 151 



u/f = i + ih + o{e) with e = 



(47) 



The value £4 — — ln(M 2 /A 2 ) depends on the intrinsic 
QCD scale A and in j^y it is suggested to use Aw 4 it f n ; 
in Ref. |§2 a value of £4 w 4.0 ± 0.6 is quoted. 

ChPT also relates the decay constant to the scalar form 
factor radius via 



1 
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U/f = l + -(r 2 ) s M^ + -^ + 0(e) 



(> 



12 



(48) 



The relation holds for both, the pion and the kaon; this 
behavior is confirmed by our results in Fig. 1131 
Eq. I|48|) may be translated to 



(r 2 )s 



(49) 



The authors of Ref. [52( quote an expected value of 
U/f = 1.072(4) and (r 2 ) s = 0.61(4) fm 2 . 

In our quenched QCD case one expects correction 
terms with a logarithmic singularity in the valence quark 
mass in. As pointed out in [4l| the leading order logarith- 
mic term m log rn of ChPT involves quark loops that are 
absent in the quenched case. There will be non-leading, 
e.g., logarithmic, terms, though. We therefore allow in 
addition to the linear term in the quark mass m also a 
term m 2 logm in the extrapolating fit (cf. the discus- 
sion in ^3)- Actually, in the fit it makes no significant 
difference whether we take this term or just m . 

When removing the leading scale dependent factors, 
thus changing to physical values, in Fig. ^] we find still 
non-negligible a-dependence away from the chiral limit. 
Note that the physical values of f„ and Jk have to be 
read off at the respective different values of mpg. 

Fig. El ( ne gl ec ti n g possible quenching effects) may be 
analyzed according to the ChPT expansion l|47|) . We see, 
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however, that the slope near the chiral limit shows con- 
siderable a-dependence. Taking the slopes at face value 
we obtain values for {r 2 ) s of 0.08 fm 2 and 0.13 fm 2 , con- 
siderably smaller than expected for full QCD. 

For a better study of the scale dependence we plot in 
Fig.^Jdata for all lattice spacings, always for the largest 
volumes studied here. This plot combines different phys- 
ical volumes. As can be seen in the upper part of the 
figure, the expected leading a-dependence is nicely ex- 
hibited. In the middle part of the figure we have removed 
the leading a-dependence. 

As mentioned above in our discussion of the conden- 
sate, we expect the linear corrections to the leading scal- 
ing behavior to be small, but cannot exclude them, since 
the action is not an exact GW-operator. Like for the con- 
densate we estimate the systematic error by fitting the 
data to a constant (no scaling violations) and a linear be- 
havior in the lattice spacing and quote the difference as 
systematic error. The constant extrapolation for the pion 
gives = 96(2) and the linear one gives /„■ = 100(10). 
We therefore obtain as final continuum extrapolation 



/ 7r = 96(2)(4) MeV 



(50) 



For the kaon we have results for non-degenerate quark 
masses only at two lattice spacings (same physical vol- 
ume) and an extrapolation is thus underdetermined. At 
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the other values of the lattice spacing studied we there- 
fore compute values for fpg from the /^-results at that 
mass to, where at the pseudoscalar state for equal-mass 
quarks agrees with the kaon mass. This is the usual 
method when only one quark mass is considered and is 
justified from the universal behavior observed, e.g., in 
Fig. 1131 The corresponding numbers are indicated by 
triangles in Fig. 1151 Extrapolating in the same manner 
as discussed for the pion we end up with 



f K = 106(1) (8) MeV . 



(51) 



We expect the leading corrections to scaling to cancel 
in the ratio Jk/J-r, plotted in the lower part of Fig. ^3 
and indeed the ratio is compatible with constant behav- 
ior. The continuum extrapolation along the lines dis- 
cussed above gives the result 



f K /U = 1.11(1)(2) . 



(52) 



Studies for the overlap action in quenched simulations 
have obtained similar results for f„ [ll HIM El BE!. 
The experimental values are / w = 92.4(0.3) MeV and 
fx = 113.0(1.3) MeV [3lJ (taking into account the factor 
v2 in the definition) . 



(at scale 2 GeV). We end up with the following final 
(physical, renormalized) values: 



Quark masses: 

Condensate: 
Decay constants: 



to = 4.1(2.4) MeV , 

to s = 101(8) MeV , 

E = -(286(4) (32) MeV) 3 

U = 96 (2) (4) MeV, 

f K = 106(1) (8) MeV , 

Sk/U = 1.11(1)(2) . 



Like many quenched results, these numbers are surpris- 
ingly close to the experimental values |31j. 
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V. CONCLUSION 



Within the Bern-Graz-Regensburg collaboration two 
Ginsparg- Wilson type Dirac operators (FP fermions and 
the so-called chirally improved Dirac operator) have been 
studied, both obeying the GWC to a good approxima- 
tion. For the chirally improved Dirac operator we know 
the quark bilinear renormalization constants |27| . This 
enables us to determine basic low energy constants for 
the light- and strange quark sector. All the results pre- 
sented here are in the quenched approximation, i.e., with- 
out taking into account dynamical sea quarks. Due to the 
complicated and expensive form of the Ginsparg- Wilson 
type operators, a full QCD study still is a task for future 
work. 

Since the chirally improved operator is substantially 
less expensive in computational means than the over- 
lap operator, it was possible to work at several lattice 
spacings ranging from 0.15 fm down to 0.08 fm and with 
different lattice sizes. The results presented here are 
based mainly on simulations for lattice spacing 0.15 fm 
and 0.12 fm on lattices with a physical spatial extent of 
2.4 fm. For an investigation of the scaling behavior and 
for studying finite size effects we also consider results 
from the other lattice parameter sets with different size 
and lattice constants. The data cover pion masses down 
to 330 MeV. 

In the range of our simulation parameters we cannot 
reliably identify behavior specific for QChPT. Therefore 
the data have been extrapolated towards the chiral limit 
with the leading order ChPT expansion for full QCD. All 
physical values have been converted to the MS-scheme 
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APPENDIX: TABLES 



Here we collect our main results for the meson and 
quark masses, for the condensate and the decay con- 
stants. All quantities are given as renormalized ones, i.e., 
including the renormalization factors converting them to 
the MS-scheme at /j, = 2 GeV. 
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TABLE IV: The pion and kaon masses computed from Eq. I|26[l on the 16 3 x 32, p = 7.90 lattice. As soon as the light quark 
mass exceed the (fixed) strange quark mass, the kaon becomes lighter than the pion. 



a m 


a 2 Ml 


Ml[GeV 2 } 


M w [MeV] 


a 2 Ml 


Mi[GeV 2 ] 


M K [MeV] 


0.02 


0.062(3) 


0.110(5) 


332(7) 


0.163(2) 


0.290(3) 


538(3) 


0.03 


0.092(2) 


0.163(4) 


404(5) 


0.178(2) 


0.317(4) 


563(3) 


0.04 


0.121(2) 


0.216(4) 


465(4) 


0.194(2) 


0.344(4) 


586(4) 


0.05 


0.151(2) 


0.269(4) 


519(4) 


0.209(2) 


0.371(4) 


609(4) 


0.06 


0.181(2) 


0.322(4) 


568(4) 


0.224(3) 


0.398(5) 


631(4) 


0.08 


0.241(3) 


0.429(5) 


655(4) 


0.254(3) 


0.452(5) 


672(4) 


0.10 


0.301(3) 


0.536(6) 


732(4) 


0.284(3) 


0.505(5) 


711(4) 


0.12 


0.362(3) 


0.644(6) 


803(4) 


0.314(3) 


0.559(6) 


747(4) 


0.16 


0.488(4) 


0.868(7) 


931(4) 


0.375(3) 


0.666(6) 


816(4) 


0.20 


0.618(4) 


1.098(7) 


1048(3) 


0.436(4) 


0.775(7) 


880(4) 




TABLE V: 


_L IlfcJ JJ1UI1 cLLlLl KdUIl 


masses computed from Eq. (12(i> on the 20 3 


y QO ft — file latfirp 




a m 


a 2 Ml 


M 2 [GeV 2 ] 


M n [MeV] 


a 2 M 2 K 


Mi[GeV 2 ] 


M K [MeV] 


0.017 


0.052(1) 


0.143(4) 


378(5) 


0.108(2) 


0.295(5) 


543(4) 


0.02 


0.060(1) 


0.164(4) 


405(4) 


0.112(2) 


0.305(5) 


553(4) 


0.025 


0.073(1) 


0.199(4) 


447(4) 


0.118(2) 


0.323(4) 


568(4) 


0.03 


0.086(1) 


0.235(4) 


484(4) 


0.125(2) 


0.341(4) 


584(4) 


0.04 


0.112(1) 


0.305(4) 


552(4) 


0.138(1) 


0.376(4) 


613(3) 


0.06 


0.164(1) 


0.448(4) 


669(3) 


0.164(1) 


0.448(4) 


669(3) 


0.08 


0.217(1) 


0.592(4) 


769(2) 


0.190(1) 


0.519(4) 


721(3) 


0.10 


0.271(1) 


0.739(4) 


860(2) 


0.217(1) 


0.592(4) 


769(2) 


0.12 


0.326(1) 


0.890(4) 


943(2) 


0.243(1) 


0.664(4) 


815(2) 


0.16 


0.439(2) 


1.200(5) 


1096(2) 


0.297(1) 


0.810(4) 


900(2) 



TABLE VI: The light quark condensate as derived from different observables (lattice size 16 3 x 32, p = 7.90). In the last line 
we give the extrapolation to the chiral limit as discussed in the text. 



a m 


a 3 |S(a 4 p) I 


a 3 ^GMOR 


a 3 |£(pp)(a 4 a 4 >| 


^/|E (A4P) |[MeVj 


{/|S GMO R|[MeV] 


^/|E (PP)(A4A4) |[MeVj 


0.02 


0.0101(5) 


0.0101(8) 


0.0104(5) 


289(5) 


288(8) 


291(4) 


0.03 


0.0102(5) 


0.0097(7) 


0.0102(5) 


289(4) 


284(7) 


290(5) 


0.04 


0.0103(4) 


0.0098(6) 


0.0103(5) 


290(4) 


285(6) 


290(4) 


0.05 


0.0105(4) 


0.0101(6) 


0.0104(5) 


292(3) 


288(6) 


291(4) 


0.06 


0.0108(4) 


0.0104(6) 


0.0107(4) 


295(3) 


291(5) 


293(4) 


0.08 


0.0113(4) 


0.0111(6) 


0.0113(4) 


300(3) 


297(5) 


299(3) 


0.10 


0.0119(4) 


0.0117(5) 


0.0119(4) 


305(3) 


303(4) 


304(3) 


0.12 


0.0125(4) 


0.0123(5) 


0.0125(4) 


309(3) 


308(4) 


309(3) 


0.16 


0.0136(4) 


0.0134(5) 


0.0136(4) 


318(3) 


317(4) 


319(3) 


0.20 


0.0146(4) 


0.0144(5) 


0.0146(4) 


326(3) 


324(4) 


326(3) 


chir. 


0.0092(5) 


0.0087(8) 


0.0092(6) 


279(5) 


274(8) 


279(6) 



TABLE VII: The light quark condensate as derived from different observables (lattice size 20 3 x 32, p = 8.15). 



a m 


a 3 |S(a 4 p>| 


O 3 pGMOR 


a 3 |E <PP ) < A 4 A 4 >| 


V\E (AiP) \\MeV\ 


^/|E GM OR|[MeV] 


^/|E < pp )( A 4 A 4) |[MeV] 


0.017 


0.0058(4) 


0.0061(3) 


0.0059(2) 


297(6) 


302(5) 


299(4) 


0.02 


0.0059(3) 


0.0061(3) 


0.0059(2) 


298(6) 


302(5) 


299(4) 


0.025 


0.0060(3) 


0.0061(3) 


0.0060(2) 


300(5) 


302(5) 


301(3) 


0.03 


0.0061(3) 


0.0062(3) 


0.0061(2) 


302(4) 


304(4) 


302(3) 


0.04 


0.0063(2) 


0.0065(2) 


0.0063(2) 


305(3) 


308(4) 


306(3) 


0.06 


0.0069(2) 


0.0070(2) 


0.0068(1) 


314(3) 


317(3) 


314(2) 


0.08 


0.0074(2) 


0.0076(1) 


0.0074(1) 


323(2) 


324(2) 


322(2) 


0.10 


0.0080(1) 


0.0081(1) 


0.0079(1) 


330(2) 


331(2) 


329(1) 


0.12 


0.0085(1) 


0.0085(1) 


0.0084(1) 


337(2) 


338(2) 


336(1) 


0.16 


0.0094(1) 


0.0094(2) 


0.0093(1) 


349(2) 


349(2) 


348(1) 


chir. 


0.0052(3) 


0.0055(3) 


0.0053(2) 


287(6) 


291(6) 


289(4) 
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TABLE VIII: Pion and kaon decay constants (from Eq. J2D) for 16 3 x 32, (3 = 7.90 lattice. In the last line we give the 
extrapolation to the (semi-)chiral limit (where the light quark masses vanish) as it is discussed in the text. 



a m 


afn/ZA 


a fir 


.MMeV] 


a jk/Za 


a JK 


jk [JvLev] 


0.02 


0.076(2) 


0.076(2) 


101(3) 


0.080(2) 


0.079(2) 


106(2) 


U.Uo 


0.076(2) 


0.075(2) 


101(3) 


U.UoU(z ) 


U.U fy(z ) 


lUO(z ) 




0.077(2) 


0.076(2) 


102(3) 




U.Uou^Z j 


1 cir(o\ 

lUO^Z ) 


0.05 


0.079(2) 


0.078(2) 


104(3) 


0.081(2) 


0.080(2) 


107(2) 


0.06 


0.080(2) 


0.079(2) 


105(2) 


0.082(2) 


0.081(2) 


108(2) 


0.08 


0.082(2) 


0.082(2) 


109(2) 


0.083(2) 


0.082(2) 


109(2) 


0.10 


0.084(1) 


0.084(1) 


112(2) 


0.084(1) 


0.083(1) 


111(2) 


0.12 


0.086(1) 


0.086(1) 


114(2) 


0.085(1) 


0.084(1) 


112(2) 


0.16 


0.089(1) 


0.089(1) 


118(2) 


0.086(1) 


0.085(1) 


114(2) 


0.20 


0.092(1) 


0.091(1) 


121(2) 


0.087(1) 


0.087(1) 


115(2) 


(semi-)chir. 


0.072(3) 


0.071(8) 


95(4) 


0.078(2) 


0.077(2) 


103(3) 




TABLE IX: Pion and kaon decay constants (from 


Eq. (HU) for 20 3 


x 32, (3 = 8.15 lattice. 




am 




a U 


WMeV] 


a /k/Za 


af K 


/x[MeV] 


0.017 


0.062(1) 


0.061(1) 


101(2) 


0.067(1) 


0.066(1) 


109(2) 


0.02 


0.063(1) 


0.062(1) 


103(2) 


0.067(1) 


0.066(1) 


110(2) 


0.025 


0.064(1) 


0.063(1) 


104(2) 


0.067(1) 


0.067(1) 


110(2) 


0.03 


0.065(1) 


0.064(1) 


106(2) 


0.0678(9) 


0.0670(9) 


111(1) 


0.04 


0.067(1) 


0.0661(9) 


109(2) 


0.0686(8) 


0.0679(8) 


112(1) 


0.06 


0.0703(8) 


0.0695(8) 


115(1) 


0.0703(8) 


0.0695(8) 


115(1) 


0.08 


0.0730(7) 


0.0721(7) 


119(1) 


0.0716(7) 


0.0708(7) 


117(1) 


0.10 


0.0752(6) 


0.0744(6) 


123(1) 


0.0728(7) 


0.0720(7) 


119(1) 


0.12 


0.0772(6) 


0.0763(6) 


126(1) 


0.0737(7) 


0.0729(7) 


121(1) 


0.16 


0.0801(6) 


0.0792(6) 


131(1) 


0.0752(6) 


0.0743(6) 


123(1) 


(semi-)chir. 


0.058(2) 


0.058(2) 


95(3) 


0.065(1) 


0.064(1) 


106(2) 
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